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<3\ ' Abstract 



Gravitational lensing is predicted by general relativity and is found in observations. When 
a gravitating body is surrounded by a plasma, the lensing angle depends on a frequency of the 
electromagnetic wave due to refraction properties, and the dispersion properties of the light prop- 
agation in plasma. The last effect leads to dependence, even in the uniform plasma, of the lensing 
angle on the frequency, what resembles the properties of the refractive prism spectrometer. The 
strongest action of this spectrometer is for the frequencies slightly exceeding the plasma frequency, 
what corresponds to very long radiowaves. 



<N ■ 1 Introduction 

_> , 

An ordinary theory of the gravitational lensing is developed for the light propagation in the 
vacuum. Gravitational lensing in the vacuum is achromatic because the deflection angle for the 
photon does not depend on the frequency of the photon pQ. In the limit of a weak lensing in 
vacuum by a body with a mass M the deflection angle for the photon (Einstein angle) is a — 
AGM/c 2 b = 2r g /b, under condition b^> r g , where b is impact parameter, r g is the Schwarzschild 
OO ! radius pQ. 

Propagation of the light in the medium at presence of the gravity field was considered by many 
authors [2], [3], [4], [5] and references there. If we consider inhomogeneous medium (without grav- 
ity) the light rays move along the curved trajectories in this medium. In the papers concerning 
gravitational deflection the inhomogeneous medium was considered. The deflection due to the 
gravitation, and the deflection due to the inhomogenity of the medium had been considered sepa- 
rately, without an account of the influence of the dispersion in plasma on the light propagation in 
the gravitational field. In this work we show that even in the homogeneous medium the dispersion 
in the plasma leads to dependence of the light deflection angle on the wavelength, what is different 
from the constant deflection angle in the vacuum. 

It have been shown [T], [5J, [B], [7], [5], that light propagation in the gravitational field in 
the vacuum may be formulated as its propagation in a inhomogeneous medium with an effective 
refraction index n gi depending on metric. It have been shown also, that in presence of the medium 
in the gravitation field, such analogy is valid too. In this case we should use the refraction index, 
which is a multiplication of the effective gravitational refraction index n g and usual refraction 
index n, determined by the physical properties of the medium: n e jf = n g n [SJ, [3]. In the case 
when both n g and n are close to unity, the combined deviation of the refraction index from unity 
is reduced to the sum of both separate effects [5] , [3] , [SJ . 

In this work we consider the gravitational lensing in a homogeneous plasma. Plasma is a 
dispersive medium, where the refraction index depends on the frequency of the photon. Therefore 
in the plasma the photons with different frequencies move with different velocities, namely the 
photons with smaller frequency (or bigger wavelength) move with smaller group velocity of the 
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light signal. We obtain here, that in a homogeneous plasma, in the presence of gravity, the 
deflection angle of the photon depends on the frequency of the photon, and discuss observational 
effects of this phenomenon. 

In the works, where the effect of the light dispersion in the plasma was not taken into account, 
the dependence of the gravitational deflection on frequency was connected only with the plasma 
inhomogeneity, and disappeared in the uniform plasma. In the book of Synge [3] , the geometrical 
optics in the medium with gravity was considered in great details, and he had derived equations 
for the photon propagation in an arbitrary medium with gravity. Here we calculate, on the basis of 
equations of Synge [9], the deflection angle for the photon wave packet, moving in the gravitational 
field in the presence of a uniformly distributed plasma. We obtain the dependence of the deflection 
angle on the frequency for the case of a week field of the Schwarzschild black hole metric. The 
deflection angle increases with decreasing of the frequency, and when the frequency is approaching 
the electron plasma frequency ui 2 = " e , all photons are falling into the black hole, if their 
impact parameter is less than the critical one, depending on frequency, b < b c (uj). When oj 
approaching u> e , the critical b c formally goes to oo. 

2 Light propagation in an inhomogeneous plasma in a week 
gravitational field 

Let us consider a static space-time with the metric 

ds 2 = g ik dx i dx k = g aSi dx a dx S3 + g 00 (dx ) 2 , i, k = 0, 1, 2, 3, a, = 1, 2, 3. (1) 
Here gik do not depend on the time. Let us assume that the gravitational field is week, so that 

9ik = Vik + hik, h ik < 1, hik — > under x a — > oo . (2) 
Here rjik is the metric of a flat space (—1, 1, 1, 1), and is a small perturbation. Note [7], that 

g ik =V ik ~h ik , v ik =Vik, h lk = h lk . (3) 

Let us consider in this gravitational field a static inhomogeneous plasma with the refraction index 
n, which depends on the space location x a and the frequency of the photon oj(x a ): 

[uj{x a )\' ! m 

Here u>(x a ) is the frequency of the photon, which depends on space coordinates x , x 2 , x 3 due 
to presence of gravitational field (gravitational red shift). We denote w(oo) = us, and e is the 
charge of the electron, m is the mass of the electron, N(x a ) is the electron concentration in the 
inhomogeneous plasma, uj e is the electron plasma frequency in this plasma. To consider a general 
case let us assume 

N{x a ) = N a + N 1 (x a ), N = const, iVi(oo) = 0. (5) 
Here N% is not supposed to be small compared to No, let us denote: 

Aire 2 

u 2 =uj 2 + w{, where iv 2 = K e N , uj\ = K e N\, K e = . (6) 

The gravitational optic in a medium in a curved space-time, was investigated in [5]. Consider a 
set of three-dimensional surfaces in a static four-D space-time, which are characterized by mono- 
tonically increasing phase angle. These 3-D surfaces are called 3-D waves, or phase waves. It was 
found in [5], the connection between the phase velocity u, and a 4- vector of the photon momentum 
p l , written with using the refraction index of medium n, n = c/u, c is the light velocity in a 
vacuum, as 

—z = n = 1 H o • (7) 
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The refraction index n, defined for plasma in is in a general case the function of x % and ui(x a ), 
which are determined by the properties of the medium, and the photon frequency. In the case 
of the vacuum (n = 1) we can obtain from (|7|) the usual relation for the square of the photon 
4-vector: pip 1 = 0. In the medium, the square of the photon 4-vector is not equal to zero. For the 
medium in a flat space-time we have, 

300 = "I, 9aa = 1, p" = -P0, P a = Pa 



-(p ) 2 + (p a ) 2 
(P ) 2 



n = 1 + Lcm . (8) 



and obtain the usual relation between the space and time components of the 4-vector of the photon 

ES> m. na> 

(p a ) 2 = n 2 (p ) 2 . (9) 

The time component of the photon 4-vector is its energy, therefore p° is proportional to the 
frequency of the photon u> [8] . We have in the flat space-time 

p° = Clj, C = const > 0, (10) 
and in a space-time with gravity we have [E] , [H] , 

P°V=g^=Cw(x a ), (11) 

what physically determines the gravitational red shift. The coordinate x a and the momentum p a 
are connected by the relation dx a /d\ = p Q , where A is a parameter changing along the photon 
trajectory (see below). Consider a photon moving along z-axis, with the frequency at infinity equal 
to oj. Without the gravity and medium inhomogeneity its unperturbed trajectory is a straight line 
along z axis. The photon 4-vector for the unperturbed trajectory is p l = (p°,0, 0,p 3 ), and the 
relation between the space and time components, using the flat metric rjik, is 

(P 3 f = nt(p°f- (12) 

Here we denote ((6]) 



no = »(oo) = yi-^. (13) 

It is convenient to use the coordinate z as the parameter A. Then we have the components of the 
photon 4-vector in a simple form: 

^ = (1/710,0,0,1), Pi = (-l/no, 0,0,1). (14) 
Thus the photon momentum in plasma is a time-like 4-vector 

P i Pi = i - \ = i - ( i - 4) < °> ( is ) 

While the phase velocity u — ^ in plasma is larger than the light velocity in vacuum c, the group 
velocity v gr is less than c, so that the larger frequency corresponds to the larger group velocity, 
tending to c in the limit. For the group velocity we have the relation [10j . [S] 



c 

v 



9'' 



V/1 -«§/w a ' 9r V w 



' 'i_^ =cno<Ci (16) 



so that v gr u = c 2 . To find a constant C in pip , we consider this relation at infinity, where 
u(x a ) = w, and p° is defined by (fT4|) . We have then 

p° = Cu, C = — . (17) 

TLqUJ 
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The trajectories of the photon in presence of the gravitational field may be found from the varia- 
tional principle [5] 



Pi dx l = 



with the restriction l[7|). which may be written in the form 



W{x\ Pi ) 



= 0. 



(18) 



(19) 



Here we define the scalar function W(x t ,pi) of x l and pi. The variational principle (fTH)) , with the 
restriction condition W(x l ,pi) = 0, leads to the system of differential equations [jj]: 



dx i _ dW dp, _ dW 
dX dpi ' dX dx 2 
with the parameter A changing along the light trajectory. Let us introduce a variable 



X = Po^9™ = -p"V^ = -Cu{x a ) = -—o;(x a ), 

UqUJ 

by using of which we can transform W(x l , P i) to a simpler form 



W{x\ Pl 



g^p l p 3 -(n 2 -l){poV Z ?°y 



(20) 



(21) 



(22) 



= \ [g oa Pm + g a0 p aPfJ - (n 2 - i) P 2 (-g oa )} = \ [g^ PaP , - »V] 

From (|20|) we obtain the system of equations for the space components p a : 



^ xOL a/3 dp a 1/37 ^ / 2 2\ 

For the inhomogeneous plasma n = n(x,x a ), and we have, using (JH),©, (1^1]) the relation 



2 ^ 



K e N K e Ni(x a ) 



X 



X 



;K e N Q 



[Uj{x a )] 2 [uj(x a )Y 

d X 1 K e dN^x") 



(23) 



(24) 



nZui 2 



dx a 2 n\uj 2 dx a 



As follows from (|20|) . the variable po is constant along the trajectory, so from (|21|) we have 

dx 



dx a 



Po 



and obtain finally 



~ (»V), a = XPO {V^) a - \ il, = ^ 



Using 



we reduce equations (|23p to the form 



<9Aq(a; Q ) _ 1 
2 



1 1 



2uj 2 -oj 2 



(25) 

9Ai(x Q ) 
dx* 

(26) 



< ^ a;Q a/3 1 Q~i , 1 2 i' OCA 1 



1 1 



A, 



dNi_(x a ) 



2 lo 2 ~uj 2 e dx a 



(27) 



The system (|27p describes the light propagation in an inhomogeneous plasma, with account the 
dispersion, in the presence of a gravitational field. Let us consider the photon moving along z-axis 
in a flat space with a homogeneous plasma, and use the coordinate z as the parameter A. In 
approximation of small perturbation hik and week inhomogeneity, Ni <C Ao, one can integrate 
equations , calculating right-hand side of equations (f2"7) by using the unperturbed trajectory of 
the photon, with p l from HH), in the right hand side. At this simplification, the first and second 
terms in the second equation (1137(1 can be written as 
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1 H%P0i + \pI (-* 00 ) « = -I^PI + \vl {-9 m ), a ■ (28) 



In the week gravitational field we can write, using |2J),([3|) and (11411 



(29) 



Finally we obtain the equation describing the light propagation in a weekly inhomogeneous plasma, 
in the presence of a small gravitational field, as 

d Pa 1 / Lo 2 h 00 \ 1 K e dNijx") 

- o ^33 + -o -o - o 7~i2 "3 ' (• 3U ) 



2 \ a; 2 — lu 2 ) 2 w 2 — lu 2 9x c 

The deflection angle in the plane, perpendicular to the unperturbed light trajectory, is determined 

as 



dtp = [P/9 (+oo) -pp{-aa)]/p, p = ^Jp\ + p% + p% = |p 3 | = 1, /3=1,2. (31) 

Here p is defined by the unperturbed trajectory, and (3 = 1, 2 are related to x, y axes. After 
integration, we have the following expression for the deflection angle of the photon with the 
unperturbed trajectory along z axis 

oo 

1 f d ( h m uj 2 K e Nx{x a ) 



&t) = 2 / ' k33 + 772 - ' 2 ~ _ 2 — ~ ' dz ' (32) 



If the problem is axially symmetric, it is convenient to introduce the impact parameter b. Let us 
consider the photon moving along the axis z in the plasma with the impact parameter b relative to 
the point mass, and the plasma has a spherically-symmetric distribution of a concentration around 
this point mass, N = N(r) = No + N\(r). The situation is axially symmetric, so the position 
of the photon is given by b and z, and the absolute value of the 3-radius- vector is r = \Jb 2 + z 2 , 
instead of r = \J x\ + x\ + z 2 . We have the following expression for the deflection angle in the 
plane perpendicular to direction of the unperturbed photon: 

oo 

1 f fdh 33 , Lo 2 dh 00 K e dNx^) \ , .... 

ab = 2j [-W + ^4^b~ u 2 - u 2 db ) dz - (33) 

— oo 

Note that db < corresponds to bending of the light trajectory to the direction of the gravitation 
center, and on, > corresponds to the opposite deflection. 

3 Particular cases 

3.1 Vacuum, and a homogeneous medium without dispersion 

Let us consider a homogeneous medium without dispersion, with the refraction index ri=const > 1 
not depending on the frequency to. The equations describing the light propagation follow from 
CZ3): 

d% a uB dp a 1 «« n 2 , 2 . 

-a =g p * -d\=-2 9 >^ + Y( x l<* (34) 

The right-hand side of equation (j2"3"|) for the p a is transformed, using (|2"Tj) , (fT4"]) , and for deter- 
mination of the deflection angle ctb of the photon, moving with the impact parameter b in the 
spherically symmetric gravitational field, we have 

-\&PfOh + y (x 2 ), tt = -\f&PfOh + yPo = ~\a>l + y p3 (-.9 00 ), Q = (35) 
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= g { h 33), a p 3 + Y Po ( h °°),c, = 5 (^aj.c + y -J ( h oo) )Q = 2 ^ ("oo + /133) 

After integration we obtain the expression for the deflection angle of the photon which coincides 
with the well known formula [13], [13], P~5] for the vacuum case, derived from the geodesic equation, 

oo 

1 f d 

otp = 2 / ^3 (^33 + h 00 ) dz. (36) 

— oo 

We see here, that in case of the non-dispersive medium, the constant index of the refraction is 
canceled, and the photon trajectory is the same as in the vacuum, in presence of the gravitational 
field, in spite of lower velocity of the light propagation in the medium. Note that the motion of 
photons in 4-space, in the medium, is not described by the geodesic equation (neither massive not 
zero), because the light propagation in the medium with the refraction is determined not only by 
the gravitational field, but also by the medium. 

3.2 Homogeneous plasma in a week gravitational field 

Here we represent the main new result of this work: the dependence of the deflection angle 
on the photon frequency in the homogeneous plasma, due to account of the dispersion. In the 
homogeneous plasma from ([53")) , we have expression for the deflection angle as 

d_ f 1_ 

db [ 33 + 1 - u*/u> 2 



& b= I si ( ^33 + q -27-^00 ) dz. (37) 





Let us calculate the deflection angle for the photon moving in the homogeneous plasma in the 
Schwarzschild metric [7J of the point mass M 

rlr 2 

ds 2 = -(1 - rjr) dt 2 + — + r 2 (d8 2 + sin 2 9 dip 2 ), (38) 

1 — r g /r 

where r g — 2G c ^ 1 is the Schwarzschild gravitational radius. In the week field approximation this 
metric is written as [7J 

ds 2 = ds 2 Q + r -?-(c 2 dt 2 + dr 2 ), (39) 
r 

where ds^ is the flat part of metric {ds^ = r]ikdx l dx k ). The components hik are written in the 
Cartesian frame as [7J 



h o = — , h a/3 = —n a n Pl h 33 = cos 2 9. (40) 
r r r 

Here n a is the unit vector of 3-radius- vector r a = (xi,X2,Xs), the angle 9 is the polar angle 
between 3-vector r a and z-axis, and cos 9 = z/r = zj\Jb 2 + z 2 . The integrals in (J37J are taken 
analytically: 

i hm dz = I ivWTW dz = ~ r » I WTzW 2 dz = ~ r » WWTT 2 + const ' (41) 

2 r q-A ^3 



d , f d r„ z 2 f 3z b 



z 



db h 33 dz - J m b2 + Z 2 dz - T 9 J {b 2 + z 2 )5 /2 r 9 b (& 2 + ^2)3/2 + C0IiSL 

(42) 

Substituting the limits of integration, we obtain the following expression for the deflection angle 
in the homogeneous plasma 

oo 

&b = I i p 3 + i-ln^) dz = - r i{ 1 + 1-li/u 2 ) ■ (43) 
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Here < for lu > u>q, what means that the light ray is bent to the direction of gravitation 
center, as it occurs in the vacuum. The formula (143|) is valid only for lu > luq, because the waves 
with lu < loq do not propagate in the plasma [12j . In the theory of the gravitational lcnsing the 
deflection angle is usually defined as the difference between the initial and the final ray directions 
a = ei n — e ou t, where e is the unit tangent vector of a ray pQ. Therefore, if we use this definition, 
we will have the expression with the opposite sign: 

which turns into the deflection angle pQ for vacuum 2r g /b, when u — > oo. For lower frequencies 
the deflection angle may be much larger than in the vacuum, and the image of the point source 
will be represented by the line, or the ring (see Fig. 1,2), on which the frequency is decreasing with 
increasing of the distance from the source on the plane of the view. Such effect may happen only 
for photons with radio frequencies, because optical frequencies are much higher than the plasma 
frequency ujq, so the effect should be negligible. We see therefore, that the gravitational lens in 
plasma is acting as the gravitational radiospectrometer, see Fig. 1,2. 



3.3 Inhomogeneous plasma with gravity 

Let us consider inhomogeneous plasma in presence of a week gravitational field. Let us consider a 
case when plasma concentration N(oo) — 0. In our notations we have Nq = 0, luq = 0, N = N±. In 
this case we have from ([55)1 the following expression for the deflection angle, valid when \n— 1| <C 1, 

oo 

. 1 f fdh 33 dh 00 1 „ dNi(r) \ 

ab =2 J \-db + ^b--* Ke ^^) dz - (45) 

— oo 

The case of the inhomogeneous plasma with gravity was considered in [2] , [3] , [5] • Let us calculate 
the deflection angle for the photon, moving in the Schwarzschild metric, with the concentration of 
plasma in the form [5] 



r 



N{r) = N m — , N m = const, R = const, h > 0. (46) 



It is evident from (|45[) that in the week gravitational field, in presence of plasma which implies 
a small perturbation in the photon propagation (lu 2 ^> lu 2 ), both effects: the deflection due 
to gravity and the deflection due to inhomogencity of medium (non-relativistic effect) can be 
considered separately. The first and the second terms in (|45|) lead to usual Einstein angle of 
deflection, following from (|4"4")) at lu 2 , =0. Let us calculate deflection due to the inhomogeneity of 
plasma, in the third term in (|45|) for the density distribution (|46|) : 



1 f ( 1 ,8#)^,_ 1 v Ar uh [d (1 



2 1 v -js^m r = -^ K ' NmR 0b [ - ]dz - (47) 



o 



After differentiation 



d ( 1 \ d 1 Kb 



db\r h J db(z 2 + b 2 ) h / 2 (z 2 + b 2 ) h / 2 + l 
we perform the integration using [16j , and properties of the T- function: 



(48) 



dz i v^r (| 



i) 



(Z 2 + fc2)V2+l hlh+\ p , 



2 / 



, T(x) = / t^e-'dt. (49) 



So we obtain 



db{z 2 + b 2 ) h / 2 b h r(|) 
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and for the third term in (14511 we have 



— oo 

Finally, for the deflection angle near the Schwarzschild metric with plasma concentration (|46|) we 
obtain the expression 

which is the same, as in [5], where this formula is written in terms of the wave length instead of 
the frequency. 



4 Discussion 

Let us consider the case of a weekly inhomogeneous medium without gravity, with the refraction 
index n = uq + n\, hq =const, n% <C uq. The system of equation for the trajectory of the photon 
in this case follows from equations (|23p with the flat metric — rjn-i see |10j . 

dx a dp a 1 g 1 , 2 2 . 1 dn 2 1 dn 

■w=f>-*r = -2 + 2 ( n * % = 2^ ^- a^- (53) 

The light propagation in a week gravitational field may be considered, as a propagation in a flat 
space with the "gravitational" refraction index, which for the Schwarzschild metric is written as 

[II HI El HI El E] 

n g = l + ^, (54) 
r r 

The total effective refraction index n e //, in presence of plasma with the proper refraction index 
n, is determined [I], [5] as n e ff = nn g . When plasma in a week gravitational field has a refraction 
index close to unity, it follows from the definition of n e f t , that in the linear approximation the 
effective refraction index is obtained as a sum of two different additions to the unity [2J, [3], [5]: 

, , r 9 ^l(r) r g uj 2 (r) 
n eff = H ^-2~, — < 1, — — < !• (55) 

From (|53|) we obtain the deflection angle (|52|) . Note that when \n 2 — 1| is not small, n e ff cannot 
be represented in the form (|55p . and such case was considered in the subsection 3.2. 

The main new result of this work is obtaining of the dependence, of the lensing angle on the 
frequency in a homogeneous plasma in the gravitational field. This effect has a relativistic nature, 
and is connected with the dispersive properties of plasma. It is interesting, that, as shown in 
the subsection 3.1, in the medium without dispersion, the trajectories of photons with different 
frequencies (energies) are exactly the same as in the vacuum, while their velocities are less that 
the vacuum light velocity c. 

Observational effect of such frequency dependence is easy to explain on the example of lensing 
by the Schwarzschild point-mass lens. This lens gives two images of source, on the opposite sides 
from lens. Angular positions of images depend on the Schwarzschild radius of lens and positions 
of source, lens and observer. The dependence of the deflection angle on the frequency in plasma 
lead to the phenomenon, that instead of two concentrated images with complicated spectra, we 
will have two line images, formed by the photons with different frequencies, which are deflected 
by different angles (Fig. 1,2). 

The description of the mass distribution as a point mass (Schwarzschild lens) is rarely sufficient 
for gravitational lensing considerations PQ . In reality the gravitational lenses have more a compli- 
cated structure, and position of images different from that of the point- mass lens. But all standard 
models of gravitational lenses are based on the same Einstein deflection angle 2r g /b, which should 
be modified for sufficiently long waves, according to our formula (|44[) . in the presence of plasma. 
Note also, that taking into account of plasma effects on the gravitational lensing may influence the 
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Figure 1: Lensing of the point source by the Schwarzschild point-mass lens. Instead of two point 
images due to lensing in the vacuum we have two line images. The pairs of images, corresponding 
to the same photon frequency, are indicated by the same numbers. Two images with number 1 
correspond to the vacuum lensing. 

spectrum of the microwave background radiation, leading to the dependence of power spectrum 
of the fluctuations on the photon wavelength. 

The light signals are propagated with the group velocity. It follows from (|4~4")) , that the smaller 
group velocity (smaller frequency and bigger wavelength) corresponds to a larger deflection an- 
gle. Hence the effect of difference in the gravitational deflection angles is significant for larger 
wavelengths, when lu is approaching uj e , what is possible only for the radio waves. Therefore, the 
gravitating center in plasma is acting as a radiospectrometer. The longest radiowaves are regis- 
tered [T7j in the band ~ 3 • 10 3 cm, corresponding to v ~ I0 7 Hz = 10 MHz and ui = 2ivv ~ 6 • 10 7 
sec -1 . The spectroscopic effects of lensing will be important when N e > 3 ■ 10 5 cm -3 , correspond- 
ing to 10 % difference in the lensing angles. Such electron densities may be expected around the 
supermassive black holes, or during lensing at earlier stages of the universe expansion at z > 10 3 . 
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Figure 2: Axis on lensing by the Schwarzschild point-mass lens. The case of the Einstein ring. Instead 
of a thin ring corresponding to the vacuum lensing (the inner circle of the ring) we have a thick ring, 
formed by the photons of different frequencies. 
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